An efficient edge based data structure for the compressible Reynolds-averaged Navier–Stokes equations on hybrid unstructured meshes

An efficient edge based data structure has been developed in order to implement an unstructured vertex based finite volume algorithm for the Reynolds-averaged Navier–Stokes equations on hybrid meshes. In the present approach, the data structure is tailored to meet the requirements of the vertex based algorithm by considering data access patterns and cache efficiency. The required data are packed and allocated in a way that they are close to each other in the physical memory. Therefore, the proposed data structure increases cache performance and improves computation time. As a result, the explicit flow solver indicates a significant speed up compared to other open-source solvers in terms of CPU time. A fully implicit version has also been implemented based on the PETSc library in order to improve the robustness of the algorithm. The resulting algebraic equations due to the compressible Navier–Stokes and the one equation Spalart–Allmaras turbulence equations are solved in a monolithic manner using the restricted additive Schwarz preconditioner combined with the FGMRES Krylov subspace algorithm. In order to further improve the computational accuracy, the multiscale metric based anisotropic mesh refinement library PyAMG is used for mesh adaptation. The numerical algorithm is validated for the classical benchmark problems such as the transonic turbulent flow around a supercritical RAE2822 airfoil and DLR-F6 wing-body-nacelle-pylon configuration. The efficiency of the data structure is demonstrated by achieving up to an order of magnitude speed up in CPU times.     

Decay rate for viscoelastic wave equation of variable coefficients with delay and dynamic boundary conditions

In this paper, we consider a viscoelastic wave equation of variable coefficients in the presence of past history with nonlinear damping and delay in the internal feedback and dynamic boundary conditions. Under suitable assumptions, we establish an explicit and general decay rate result without imposing restrictive assumption on the behavior of the relaxation function at infinity by Riemannian geometry method and Lyapunov functional method.     

改进型单晶炉提拉系统偏心平衡研究

随着光伏行业的快速发展, 对硅单晶的品质和长晶装备的稳定性的要求也不断提高。直拉法是生产硅单晶的主要方法,通过提高单晶炉副室的高度以扩大单晶硅的生产规模。由于副室高度的大幅增加,且单晶炉提拉头质心相对于旋转轴心有一定距离,对单晶炉整体稳定性有较大影响,从而降低了单晶硅的生产质量。针对此问题,对单晶炉建立可靠的力学分析模型,采用数值仿真方法,对单晶炉整体进行动力学响应分析,计算得到副室高度增加后的单晶炉工作时中钨丝绳下端晶棒的运动规律以及最大摆动幅度,为改进设计提供依据。数值仿真分析表明提高单晶炉副室高度后,提拉头较大的质心偏心是单晶炉提拉系统发生摆动的主要原因。在此基础上提出在提拉头上添加质心调节装置,通过控制系统调节可保证提拉头质心位置在旋转轴线上以降低提拉系统的摆动。     

梁柱简化模型在自立塔塔线体系动态响应中的适用性研究

对输电塔进行合理简化可以提高塔线体系动力学仿真的效率。本文给出自立塔梁柱简化模型的计算方法，并提出利用梁柱简化模型计算方法建立自立塔塔线体系整体模型，同时采用桁梁混合模型建立精细化塔线体系整体模型，对两种模型塔线体系静力特性及振型和固有频率等动力特性进行对比分析。以脱冰工况为例，采用生死单元技术将施加在输电线节点上的集中质量单元杀死来模拟脱冰，实现对塔线体系动力学响应的有限元模拟，研究塔线体系简化模型在动态响应中的适用性。结果表明，两种模型弯曲变形误差小，低阶的振型相同，固有频率值误差小，动力特性基本相同；脱冰工况下，自立塔节点位移和塔材内力时程曲线一致，在提高计算效率的情况下，能有效保证计算精度。     

Versatile anisotropic mesh adaptation methodology applied to pure quantity of interest error estimator. Steady,laminar incompressible flow

We introduce a new flexible mesh adaptation approach to efficiently compute a quantity of interest by the finite element method. Efficiently, we mean that the method provides an evaluation of that quantity up to a predetermined accuracy at a lower computational cost than other classical methods. The central pillar of the method is our scalar error estimator based on sensitivities of the quantity of interest to the residuals. These sensitivities result from the computation of a continuous adjoint problem. The mesh adaptation strategy can drive anisotropic mesh adaptation from a general scalar error contribution of each element. The full potential of our error estimator is then reached. The proposed method is validated by evaluating the lift, the drag, and the hydraulic losses on a 2D benchmark case: the flow around a cylinder at a Reynolds number of 20.     

A control volume finite element method for three-dimensional three-phase flows

A novel control volume finite element method with adaptive anisotropic unstructured meshes is presented for three-dimensional three-phase flows with interfacial tension. The numerical framework consists of a mixed control volume and finite element formulation with a new P₁DG-P₂ elements (linear discontinuous velocity between elements and quadratic continuous pressure between elements). A “volume of fluid” type method is used for the interface capturing, which is based on compressive control volume advection and second-order finite element methods. A force-balanced continuum surface force model is employed for the interfacial tension on unstructured meshes. The interfacial tension coefficient decomposition method is also used to deal with interfacial tension pairings between different phases. Numerical examples of benchmark tests and the dynamics of three-dimensional three-phase rising bubble, and droplet impact are presented. The results are compared with the analytical solutions and previously published experimental data, demonstrating the capability of the present method.     

An improved depth-averaged nonhydrostatic shallow water model with quadratic pressure approximation

Phase-resolved information is necessary for many coastal wave problems, for example, for the wave conditions in the vicinity of harbor structures. Two-dimensional (2D) depth-averaging shallow water models are commonly used to obtain a phase-resolved solution near the coast. These models are in general more computationally effective compared with computational fluid dynamics software and will be even more capable if equipped with a parallelized code. In the current article, a 2D wave model solving the depth-averaged continuity equation and the Euler equations is implemented in the open-source hydrodynamic code REEF3D. The model is based on a nonhydrostatic extension and a quadratic vertical pressure profile assumption, which provides a better approximation of the frequency dispersion. It is the first model of its kind to employ high-order discretization schemes and to be fully parallelized following the domain decomposition strategy. Wave generation and absorption are achieved with a relaxation method. The simulations of nonlinear long wave propagations and transformations over nonconstant bathymetries are presented. The results are compared with benchmark wave propagation cases. A large-scale wave propagation simulation over realistic irregular topography is shown to demonstrate the model's capability of solving operational large-scale problems.     

Interfacial Synthesis of a Monolayered Fluorescent Two-Dimensional Polymer through Dynamic Imine Chemistry

A fluorescent monolayered two-dimensional polymer (2DP) containing both tetraphenylethylene (TPE) and imine linkages is synthesized at air-water interface using the Langmuir-Blodgett method. We designed TPE-based monomers with long distances between the TPE and the imine linkages to avoid the charge transfer and therefore keep the fluorescence. A monolayered 2DP provided with more than 104 μm² in domain size and around 0.8 nm thickness was obtained through a successive Schiff base reaction at air-water interface. The nanostructures and fluorescent property of 2DP films were characterized by optical microscopy, SEM, TEM, AFM and fluorescence spectrum. Most importantly, the tip-enhanced Raman spectroscopy (TERS) was utilized here to confirm the success of the polycondensation of monolayered 2DP.     

基于MRS-SJC-Copula模型对A股与港股的动态联动性研究

由于经济增长周期变化,导致不同股票市场存在高低状态转换的现象,在研究不同股票市场之间联动性的研究时,需要考虑股票波动均值和方差两种结构变化。基于具有马尔可夫状态转换的动态SJC-Copula,结合修正ICSS算法对我国内地股票市场和香港股票市场之间的联动性进行方差结构突变点的检验。实证结果表明:国内和香港股票市场之间存在非线性非对称的时变相依性,并持续存在高低两种不同状态的概率转换。股票指数由于动态联动受到负面消息的下跌幅度大于正面消息的变化幅度,且上下尾部均受上期信息的持续影响。“沪港通”、“深港通”、中美贸易战等因素使得其上下尾部发生结构突变,内地与香港股票市场的联动性增大和市场波动幅度趋强。